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Abstract. In this paper we study one dimensional parabolic free bound- 
ary value problem with a nonlocal (integro-differential) condition on the 
free boundary. We establish global existence-uniqueness of classical so- 
lutions assuming that the initial-boundary data are sufficiently smooth 
and satisfy some compatibility conditions. Our approach is based on 
analysis of an equivalent system of nonlinear integral equations. 
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1. Introduction 

In this paper we consider the following free boundary value problem. 
Problem P. Find s(t) > and u(x,t) such that 

(1.1) Ut = u xx — \u, A = const > 0, < x < s(t), t > 0; 

rs(t) 

(1.2) s'(t)= / (u(x,t) - a) dx, t>0, a = const > 0; 

Jo 

(1.3) u(0,i) = /(t), f(t)>0, t>0; 

(1.4) u(x,0) = (p(x), <p(x)>0, x€[0,6], s(0) = 6>0, </>(0) = /(0); 

(1.5) u x (s(t),t) = 0, t>0. 

This is one-dimensional free boundary problem with unknown boundary 
x = s(t). Notice that (jl.2p is a nonlocal condition on the free boundary, and 
(jl.3p - (jl.5p are mixed type boundary conditions for the parabolic equation 

CD. 

The aim of this paper is to investigate the existence and uniqueness of 
classical solutions of Problem P. To this end we reduce the problem to a 
system of nonlinear integral equations and analyze its local solvability. The 
same approach has been used in many papers on the one-dimensional Stefan 
problem and its variations in order to prove existence-uniqueness results - 
see, for instance, jTJ Ch. 8] and |12] and the bibliography given there. In 
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the context of tumor models, existence-uniqueness results for free boundary 
problems similar to Problem P are obtained in [8, Theorem 3.1] and [H 
Theorem 2.1] by the same method. However, the presence of mixed type 
boundary conditions in Problem P brings to some additional difficulties, 
and as far as we know this case has not been studied yet. 
Our main result is the following. 

Theorem 1.1. Suppose that 

(1.6) f{t) G C\[0, oo)), f(t) >Ofort>0, <p{x) G C\[0, &]), 

<p(x) >0 forxe [0,6], /(0) = ¥>(0), /'(0) = <//'(0) - A^O), <f/(b)=0. 
Then there exists a unique pair of functions u(x,t) and sit) such that 

(i) u(x,i) is defined, continuous and has continuous partial derivatives 
u x , ut, u xx in the domain {(x,i) : < x < sit), t > 0}; 

(ii) sit) e 0^10,00)); 
(Hi) (E2P-(E^ hold. 

In order to prove this theorem we introduce an auxiliary free boundary 
value problem (see Problem P in Section 3) and analyze local and global 
in time solvability of the resulting pair of free boundary value problems. In 
Lemma 13.11 it is shown that every solution of the main Problem P ( (jl.ip - 
(jl.5p ) generates a solution of the auxiliary problem and vice versa. Existence 
and uniqueness of local solutions of the auxiliary problem are proved by 
deriving and studying an equivalent system of nonlinear integral equations 
(see Lemma [4. 21 and Lemma [5.ip . In Lemma f6.1l we obtain a priori estimates 
for the local solutions of the auxiliary problem by applying an appropriate 
maximum principle (see Lemma 12.11 and Lemma I2.2p to the solutions of 
the main problem. Finally, we prove existence of global solutions for both 
the main and the auxiliary problems by using the corresponding a priori 
estimates obtained in Lemma 16.11 Some of these results are announced 
without proofs in |15j . 

2. Preliminary results 

Throughout the paper we assume that the functions / and ip satisfy the 
conditions (|1.6p . 

Definition 1. We say that a pair of functions («(», t), s(t)) is a solution 
of Problem P for t G [0, T), T < oo, if 

(i) u(x, t) is defined, continuous and has continuous partial derivatives 
u x , Ut, u xx in the domain Dt = {i%,t) : < x < sit), < t < T}; 

(ii) the equation (jl.ip is satisfied for t < T; 

(iii) sWGC^O.T)); 

(iv) the conditions (jT3^ - ([P]) hold for t G [0,T). 

Lemma 2.1. (Maximum Principle) Let X = const > 0, s(t) G C 1 ([0,T]), 
and let u(x,t) be defined and continuous in the domain Dt = {(x, t) : < 
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x < < t < T}, have continuous partial derivatives u x ,u xx for < 
x < s(t), < t < T, and have a continuous partial derivative ut for < 
x < s(t), < t < T. Suppose that 

(2.1) -T^ + ^2-Aii>0 for 0<x < s{t), < t < T, 
and 

du 

(2.2) —(s(t),t)<0 forO<t<T. 
ox 

If M := max{ii(x,i) : (x,t) G Dt} > 0, then u(x,t) attains its maximum 
only on the union of the segments {(0,t) : t £ [0. T]} and {(x,0) : x £ 

Mo)]}- 

Proof. To the contrary, assume that u(x, t) attains its maximum M > at 
a point (xo,to) with xq > 0, to > 0. 

(a) Suppose xq < s(to). Obviously, u x (xq, to) = 0. On the other hand, in 
view of (12.11). we have 



< - (u(x , t ) - u(x Q , t - h)) = u t (x , t h ) < u xx (x , t h ) - Xu(x , t h ), 
h 

with < h < to and to — h < th < to- Therefore, letting h — > 0, we get 
u xx( x O:to) ^ AM > 0. But then the function of one variable u(x,to) has a 
strict local minimum at x = Xq, which is impossible. 

(b) Assume that xo = s(to). Since s'(to) exists, there is a unit vector 
I = (a, (3) with a > and j3 > such that the segment {(s(to) — ah, to — 
/3h), h E (0, e)} is in the interior of Dt for sufficiently small e > 0. Then 

! du 

< t (u(s(t ),t ) - u(s(t ) - ah,t - Ph)) = —=,(x h ,t h ), 
ft o£ 

where x h = s(t ) - a8h, t h = t - {39h, 9 = 9(h) £ (0, 1). Therefore, by 

(1270) it follows 



< au x (x h ,t h ) + f3u t (x h ,t h ) < au x (x h ,t h ) + f3(u xx (x h ,t h ) - \u(x h ,t h )). 

Passing to a limit as h — > we get 

< au x (s(to),t ) + (3(u xx (s(t ),to) - Xu(s(t ),t )). 

By (|2.2p we have u x (s(to), to) < 0. Therefore, the latter inequality yields 
u xx (s(to), to) > AM > 0, which implies that u x (s(to) — h, to) < for all 
sufficiently small h > 0, say h £ (0,S). But then it follows that u(s(to) — 
h, to) > u(s(to),to)) = M for h £ (0,(5), which is impossible. This completes 
the proof. □ 



In view of (jl.2p . the Maximum Principle yields immediately the following 
a priori estimates for u(x,t) and s(t). 

Lemma 2.2. If a pair of functions (u(x,t), s(t)) is a solution of Problem P 
for < t < T < oo, then 

(2.3) < u(x, t) <C T , 0<x< s(t), < t < T, 
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(2.4) - as{t) < s'{t) < (C T - a)s(t), be~ at < s(t) < be { ° T ^ ] \ 

where b = s(0), Ct = max ^sup [o.t) su P[o,b] fi x ) 

3. Auxiliary free boundary problem 
Next we consider the following auxiliary free boundary value problem. 
Problem P. Find s{t) > and u(x,t) such that 

(3.1) u t = u xx , 0<x<s(t), t>0; 

(3.2) s '(t) = ( f(t)-a)s(t) + e~ xt J (j u(£,t)dt) dx, t > 0; 

(3.3) u x (0,t) = f(t), /(i)eC([0,(X))); 

(3.4) u(x,0)=<p(x), <p(x) G ^([0,6]), s(0) = b, <p'(Q) = /(0), ^(6) = 0; 

(3.5) u(s(i),i) = 0, i>0. 

Definition 2. We say that a pair of functions u(x, t) and s(t) is a solution 
of Problem P for t € [0, T), T < oo, if 

(i) u(x,t) is defined and continuous in the domain = {(x,t) : < 
x < s(t), < t < T}, has continuous partial derivative tt x in Dt, and has 
continuous partial derivatives ut, u xx for < x < s(t), < t < T; 

(ii) the equation (13. ID is satisfied for t G (0,T); 

(hi) s^ecHlo.r)); 

(iv) the conditions J32J)-(|33J| hold for t G [0,T). 

The next lemma gives the relation between Problem P and Problem P. 
Lemma 3.1. Let f(t) and <p(x) satisfy tl.6\) . A = const > 0, and let 

(3.6) f(t) = j t (>/(*)) , <^W(*). 

(a) If a pair of functions u(x,t) and s(t) is a solution of Problem P for 
t G [0, T), then the pair of functions u(x,t) = e xt u x (x,t) and s(t) is a 
solution of Problem P for t G [0, T). 

(b) If a pair of functions u(x,t) and s(t) is a solution of Problem P for 
t G [0,T), then the pair of functions (u(x,t), s(t)) with 

(3.7) u(x,t) = f(t) + e~ xt f u(C,t)d( 

Jo 

is a solution of Problem P for t G [0, T). 

Proof, (a) Notice that u(x, t) is a C°°-function in the interior of the do- 
main Dt due to general smoothness theorems (see Ch.3, Thm. 11], and 
Corollary 2 there). Therefore, the function u(x,t) has continuous partial 
derivatives ut, u xx for < x < s(t), < t < T, and it satisfies the equation 
u t = u xx m that domain. 



GLOBAL EXISTENCE OF CLASSICAL SOLUTIONS 5 

Letting x — > in the equation (II. ip , we obtain 

u t (0,i) = f'(t) = u xx (0,t)-Xf(t). 

Thus, u x (0,t) = u xx (0,t)e xt = [f'(t) + Xf(t)]e xt , i.e., holds with f(t) = 
i {f(t)e xt ) . Now one can readily verify that the pair of functions u(x,t) = 
e xt u x (x,t) and s(t) is a solution of Problem P for t G [0,T). 

(b) We check first that the function t) given in (13.70 satisfies the 
equation (II. ID . By (j3.7l) . we have 

In order to find and justify a formula for ut we set 

t) = f{t) + e~ A * f u(£, i)d£, n = 1, 2, . . . . 

Then u n (x, t) — > u(x, t) as n — > oo for (a;, i) G -Dt- In view of (|3.1j) . we have 

O /'X f*X 

( Un (x,t)) = f'(t)-Xe~ xt / fi&tK + e"* / u K (C,t)de. 

Jl/n il/n 

Since fu n t)d£ = u x (x,t) — u x (l/n,t), we get as n — )• oo 

^(n„(x, i)) -> /'(t) - Ae" A * u(£, + e~ At (^(x, t) - fi^O, t)) 

uniformly on any compact subinterval of (0, T). Therefore, Ut(x,t) exists, 
and using (|3.7p and (|3.3j) we obtain 

ut(x, t) = /'(t) - X(u(x, t) - /(*)) + e- A *(^.(x, i) - /(£)). 

Since /(f) = 4 (f(t)e xt ) it follows that u(x,t) satisfies the equation (11. ip . 
Now one can easily see that the pair of functions (u(x,t), s(t)) is a solution 
of Problem P for t G [0, T). □ 

In view of Lemma 13-H Theorem 11.11 will be proved if we show that the 
following statement holds. 

Theorem 3.2. Suppose that 

(3.8) f(t) GC^aCoo)), f(t)>0 fort> 0, f(t) G C([0, oo)), 

£( a ;)eC 1 ([0,&]), /(0) = </5'(0), £(&) = 0. 
T/ien Problem P has a unique solution for < i < oo. 
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4. System of integral equations 

In this section Problem P is transformed to an equivalent problem of 
solving a system of nonlinear integral equations. We begin with some pre- 
liminaries. 

Consider the function 

(4.D Jf(«,«,r)- 5 ^ I _«p(-|^), r<t. 
We shall make use of the following elementary inequalities: 

* i r°° ,2 . 1 



rt -y poo 

(4.2) / \K x (x,t;£,T)\dT = -= I e~ z2 dz< 

Jo V 71 " J^=M 

(by performing the change of variable z = J^jL ); 



(4.3) [ b K(x,t;t,0)dt = -!= f^e- 

Jo V*J^ 

(by using the change of variable z = ); 

a 



(4.4) z a e-~ <z < 1-1 if z > 0, a > 0, 7 > 0. 

The next statement is a slight modification of Lemma 1 in Ch.8]), and 
its proof is the same. 

Lemma 4.1. I/s(t) 6 C 1 ([0,T]), #(t) G C([0,T]) and <t <T, then 

, s/ K x {x,t;s{T),r)g{r)dT = — — + / 2fa;(a(to),to;s(r),r)fli(r)dr, 

where in the limit we consider only points {x,t) with x < s{t). 

Next we derive a system of integral equations related to Problem P. Let 
N{x,t;£,r) be the Neumann function for the half-plane x > 0, i.e., 

iV(x, t; e, r) = i; £, r) + tf(-s, t; £, r). 

Suppose that the pair of functions {u{x, t), s{t)) is a solution of Problem P. 
For t > 0, we integrate the identity 

<"> |("*-f«)->* 

over the domain £ < r < t — e, <5 < £ < s(r) — <5, e = const > 0, 5 = 
const > 0, and pass to limits, first as 5 — > 0, and then as e — > 0. Since 
N^{x,t;0,r) = and 

rs(t-e) 

lim / N{x,t;£,t — e)u{£,t — e)d£ = u{x,i), 
z^oJo 
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it follows, in view of (13. 3p and (I3.4p . that 

5 

(4.6) u{x,t)=^J v (x,t), 
where 

(4.7) Ji(x,t)= [ N(x,t;s(T),T)v(T)dT with v{t) := u x (s(t),t), 

Jo 

J 2 (x,t) = - [ N(x,t;0,T)f(r)dT, J 3 (x,t) = [ N(x, t; £, 0)<p(t)d£ 
Jo Jo 

and 

(4.8) J 4 ( x ,t) = - [ N^(x,t;s(T),T)u(s(T),T)dr, 

Jo 

J 5 (x,t) = / s'(T)N(x,t;s(T),T)u(s(T),T)dT. 
Jo 

The condition (|3.5p implies J^{x,t) = 0, J§{x,t) = 0. Thus, the following 
integral representation holds: 

(4.9) u(x,t) = Ji(x,t) + J 2 (x,t) + J 3 (x,t). 

Next, in order to obtain an integral equation for v(t) = u x (s(t),t), we 
differentiate (|4.9p with respect to x and pass to a limit as x — > s(t) — in 
the resulting identity. In view of Lemma 14.11 h follows that 

lim ^{x,t)= f N x (s(t),t;s(T),T)v(T)dT + ^-. 

x->s(t)-0 OX Jq I 

It is easy to see that 

lim -*( x ,t) = - N x (s(t),t;0,T)f(T)dr. 
► s (t)-o dx y 



Now, consider the Green function for the half-plane x > 
G(x, t; e, r) = t; £, r) - K(-x, t; £, r). 
Since iV z = — G^, an integration by parts leads to 

(4.10) t) = - ^ G ? (x, t; e, 0)^(0^ = ^ G(x, t; £, 0)£'(£K 
because G(x, t; 0, 0) = and (pip) = 0. Therefore, it follows that 

lim -^-(x,t) = / G( S (t),t;e,0)^(C)de. 

a;-»s(t)-0 OX Jq 

Hence, for t > the function satisfies the integral equation 

(4.11) v(t) = 2[ N x (s(t),t;s{T),T)v{T)dT 
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-2 f N x (s(t),t;0,T)f(r)dT + 2 C G(s(t), t; £, 0)<ff(£)dZ. 
Jo Jo 



On the other hand, from (|3.2p and (|4.9p it follows 
(4.12) 

s '(t) = (f(t)e~ xt -a)s{t) + e- xt I / / iV(£, t; s(r), T)v{T)dTd£dx 



JO JO 



/ / N(£,t;0,T)f(T)dTc%dx+ / / / N(Z,t;r),Q)v(ri)dr]d£dx 
o jo Jo Jo Jo Jo , 



The system of nonlinear integral equations (|4.11j) and (14. 12h considered with 



s(i) = 5 + J * s'{r)dT is equivalent to Problem P, i.e., the following statement 



holds. 

Lemma 4.2. Problem P for t < T is equivalent to the problem of finding a 
pair of continuous functions (v(t), s'(t)) on [0, T) which satisfies for t > 
the system of nonlinear integral equations f^, 1 1\ ) and {4-12^ considered with 
s(t) = b + £s'(r)dr. 

Proof. We have already proved that if a pair (u(x,t), s(t)) is a solution of 
Problem P for t < T, then the pair of continuous functions v(t) = u x (s(t),t) 
and s'(t), t G [0,T), satisfies for t > the system of nonlinear integral 
equations (|4.1ip and (|4.12p considered with s(t) = b + f s'(T)dT. 

Conversely, suppose that a pair of continuous functions v (t) and s'(t), t G 
[0,T), satisfies for t > the system of integral equations (|4.1ip and (|4.12p . 
Set 



(4.13) u(x,t) 



Yl=i Jv(x, t) for < x < s(t), < t < T, 
<p(x) for < x < b, t = 0, 



where J u (x,t), v = 1,2,3 are given by gTf]) and s(t) = b + £s'(r)dT. 
We shall prove that the pair of functions (u(x, t), s(t)) form a solution of 
Problem P for t < T. 

First we show that the function u(x, t) is continuous in the domain Dt- 
Indeed, since the integrands in J\(x,t) and J2(x,t) are dominated by a 
multiple of (t — r) -1 / 2 , we have 

lim J u (x,t) = 0, v = 1,2, x £ [0,b\. 

So, it remains to show that Js(x,t) — > (p(xo) as (x,t) — > (xo,0), xq G [0,6]. 
Performing the change of variable z = (£ ± x)j2\ft, we obtain Jz(x,t) = 
Jg(x, i) + J|(x,t), where 

b — x b-\-x 

Jl(x,t) = — / 2 * e~ z2 ip(x+2zVi)dz, J|(x,t) = — e~ z2 (p{-x+2zyft)dz. 

V7TJ — 3L= ' Vvrj^i 

2\/t 2^t 



GLOBAL EXISTENCE OF CLASSICAL SOLUTIONS 9 

Now, for every xo G (0, b), it follows that 

J3 (x, t) — > (p{xo) and Jf(x, t) — >0 as (x, t) — > (xo 5 0) 

because — x /2y/i -00 and (b ± x)/2-v/i ->■ +00. 

The corner points (0, 0) and (6, 0) need a special consideration. If (x, i) — > 
(6,0), then Jf(x,t) -)• by the same argument. Since </5(6) = 0, it follows 
that Jg(x,t) — >■ as well, so J%{x, t) — > = 92(6) as (x,i) — > (b,0). 

In order to show that Js(x, t) — > y3(0) as (x, i) — > (0, 0) we shall prove that 
Jz{x n -,t n ) — > <p(0) for every sequence (x n ,t n ) — > (0,0). Since J^(x,t) and 
jf(x,i) are bounded, the sequence {Jz(x n ,t n )} is bounded. Therefore, it is 
enough to show that every convergent subsequence of the form { Js{x nk , t nk )} 
has a limit equal to <^(0). We may assume that x nk /2^/t nk — )■ a G [0, 00] 
(otherwise we may pass to a subsequence of (n^)). Then it follows that 

h{x nk ,t nk ) -> — / e~^(0)(fe + — / e-^(^(0)fiz = £(0). 

Thus, Js(x,t) — > if(0) as (x,t) — >• (0,0), which completes the proof of conti- 
nuity of u(x,t) in the domain L>t. 

It is easy to see that each of the integrals J u (x,t) is a C°°-function in 
the domain < x < s(t), < t < T, and satisfies the heat equation there. 
Thus, u(x, t) satisfies (13. ip as well. 

The functions v(t) and s'(t) are defined and continuous on [0, T) and 
satisfy the integral equations ([4. lip and (|4.12p for t £ (0,T). Consider the 
limit of the right-hand side of ([4. lip as t — )■ 0. It is easy to see that the first 
two integrals there converge to zero. With z = (£ ± s{t))/2y/t, the integral 
$ G(s(t),t;£,0)<p'(£)d£ is equal to 

b-s(t) b+s(t) 

—= / 2Vl e~ z2 <2>'(s(t) + 2z\Tt)dz = / ^ e~*V(-s(*) + 2z\ft)dz. 

As t — > 0, the first integral in the above expression tends to <p(b)/2, while 
the second one tends to zero. Therefore, by passing to limit as t — > in 
(|4.1ip we obtain 

(4.14) v(0) = (p'{b). 

Next we prove that u x (x,t) extends as a continuous function on Dt- In 
view of flUFD and (^ - (flTO]) . we have 

(4.15) u x (x,t) = h(x,t) + I 2 (x,t) + h{x,t), 
where = {(x,t) : < x < s(t), < f < T}, and 

(4.16) h(x,t)= [ N x {x,t;s{T),T)v(T)dT, 

Jo 

I 2 (x,t) = - f N x (x,t;0,r)f(r)dr, I 3 (x,t) = f G(x, t; £, 0)tf 
Jo Jo 
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We shall prove that 

(») u x (x, t) -> /(to) as (x, t) -> (0, t ), (x, t) G £>t, < t < T; 
(ii) Ua;(x,t) ->• v(i ) as (x,i) ->• (s(t ),to), (x,t) £ D?, < t < T; 
(m) Ua;(x, t) — >• <p'(xq) as (x,t) — )• (xo,0), (x,t) £ D^, < xo < b. 

Since the functions /(t), <p'(x), v(t) are continuous and /(0) = <£>'(0), <^'(&) = 
v(0) (see (|3.4p and (j4.14|) ) . the conditions (i) — (iii) guarantee that u x ex- 
tends continuous function on Dt- 

First we prove (i). Taking into account that N x (0, t, £, r) = and G(0, t, £, r) 

0, it is easy to see that I u (x,t) — > (0, to) for v = 1, 3 as x — > +0, t — > to £ 
(0,T). With the change of variable 2 = 2 Jt- T ^ e in tegral ^(x, t) becomes 

I2(x,t) = -J N x (x,t;0,T)f(T)dT = -^= J ^ e~ z fU- — \dz. 

lyft 

As x — )• +0, t — > to € (0, T) the latter integral tends to /(to). Thus (i) 
holds. 

Next we prove (ii). One can easily see that 

I„(x,t) ->■ I v (s(to),t ) as (x,t) ->■ (s(t ),t ), v = 2,3. 
From Lemma |4. II and (|4.16p it follows that 

lim Ii(x,t) = -v(to) + h(s(t ),t ). 

(x,t)-t(s(ta),to) 2 

Therefore, by (|4.1ip . we obtain 
1 3 

u x (x,t) -> -u(to) + 5Jij/(s(£o),*o) = w (*o) as ( s (*o),*o), 
i/=i 

1. e., (ii) holds. 

It is easy to verify (iii) for xo G (0, b). However, it is much more compli- 
cated to prove (iii) for xo = or xo = b. 

Next we show that (iii) holds for xo = 6, i.e., 

(4.17) u x (x,t) (p'(b) as (x, t) -)• (6, 0), (x,t)eDy. 

Let {(x n , t n )} be an arbitrary sequence such that (x n , t n ) — >■ (6,0), (x n ,t n ) G 
.D^. In order to prove that u x (x n ,t n ) — > <p'(b) it is enough to show that for 
every subsequence {(x nfe ,t nfe )} there is a sub-subsequence {(x nfc ,t nfc )} 
(which we denote for convenience by {(x m ,t m )}) such that 

Wa;(x m ,i m ) ->■ (p\b) as m -4 oo. 

We may assume without loss of generality (otherwise one may pass to an 
appropriate subsequence) that 

(4.18) — -j= — > a G [0, ooj as m — > oo. 

2\/tm 
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(Notice that x m < s(t m ) and 

b - s(t m ) s(0) - s(t m ) 

' = — - as m ->■ oo; 

therefore, every cluster point a in ()4.18j) is nonnegative.) 

In view of (|4.15p . in order to evaluate hrrim—^oo u x {x m , t m ) one needs to 
find lim m _ 5 . 00 I v (x m , t m ), v = 1, 2, 3. First, consider the case v = 1. We have 
h(x,t) = Ii t x(x,t) + h^(x,t), where 

1 /"* s f T ) _ £ (s-s(t)) 2 

(4.19) h A (x,t) = — / ; e ) u(r)dr, 



¥J 4(t-r) 3 / 2 
(4.20) / 1)2 (x,t) = 7=/ , x , /9 e" 4(t - T) v(r)dT. 



vWo 4(t-r) 3 / 2 

One can easily see that Ii,2(x,t) — > as (x,i) — > (6,0) because x + s(r) > 
6 > for (x, r) close to (6,0). On the other hand, 

(Xm, im) = -^11 (Xjtd ^m) ~l~ -^11 (%i ^m) "I - -^1,1 (%i tm)i 

where 



/ 1 /" s(t) - s(t m ) K-'Ml 2 . , , 

V 71 " 7o 4 (*m - T) ' 



r2 j. \ ^ /" m g (^m) X m 

VTTJO 4(t m -"r)'V 2 



g 4(tm-r) — g 4(t m -r) 
2 



Ii X (x m ,t m ) = — — -y^e *(*m-r) V (t)(1t. 

- /o 4(t m - Tfl 1 



v(r)dr, 



Since s'(i) is continuous, the Mean Value Theorem implies that \s(t m ) — 
s(t) I < const • |i m — t|. Therefore, the absolute value of the integrand 
of il i(x m ,t m ) does not exceed C '/ ' \Jt m — r, which leads to Iii(x m ,i m ) < 
CV^m — >• as m — )• oo. 

Changing the variable in J 3 -, by s ^™^~ Xm . — Zj we obtain 

r 3 , f X 1 f°° -z 2 ( + {S{t m )-X m ) 2 \ 1 

The expression in the square brackets in the integral I 2 1 can be written 
as exp [- {X 4^ir) )2 ) ( e9m(r) - l) . where 

9m(r) = -^TTT" "T^ • («(*m) + s ( r ) ~ 2x m) ~> as m->00 

4(t m - r) 

uniformly for r 6 [0, i m ]. Therefore, the same change of variable as in lf x 
shows that if i(x m ,t m ) — ^ as m — > oo. 
Hence, we obtain 

1 f°° _ 2 

(4.21) 7i(x m ,i m ) — >• <y(0) • —= I e z dz as m — > oo. 

V 71 " 7 Q 
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It is easy to see that 

(4.22) hixm, t m ) — > as m — > oo. 

Next we evaluate lim m _ ! . 00 Iz(x m , t m ). Since G(x,t;£,r) = K(x,t;^,r) — 
K(—x, t; £, r), performing the change of variable z = (£ T x)j2\ft we obtain 
that Is(x,t) = 13,1(2;, f) + l3,2(x,t), where 

b — x ^ b-\-x 

h,i = —}= I " * e~ z2 <p'(x+2zVi)dz, I 3 o = —= \ * * e~ z ' V ' (-x+2zVi)dz. 
V 71 " J — 2= V 71 J^= 

Therefore, in view of (|4.18|) . it follows that 

1 f° 2 

h,i(xm,t m ) -> <p'{b) ■ \ e~ z dz, h,2{x m ,t m ) ->• as m->oo, 

V 71 " J-oc 

which yields 

1 f a _ 2 

(4.23) l3(x m ,t m ) -¥ (p (b) • —j= f e z dz as m — > oo. 

V 71 " V_oo 

Hence, (HTTP follows from (l4T2TD - (|4T23j ) and IjOIjl . 
A similar argument proves that 

n 2 (x,t) -> ^'(0) as (x,t) -> (0,0), (x,t) G 

which completes the proof of (in). 

In order to complete the proof of Lemma 14.21 it remains to show that the 
condition (|3.5p holds, i.e., u(s(t),t) = for t G [0, T). Since u(x,t) satisfies 
(|3,ip - (|3.4p (as we proved above), by integrating the identity (|4.5p over the 
domain e < t < t — e, 5 < £ < s(r) — 5, s, S > 0, and passing to limits, 
first as (5 —7- and then as e — > 0, we obtain the integral representation (|4.6p . 
Now, in view of (|%B|) and (|%I5j) . it follows that 

iV € (x,t;s( 7 -),r) 5 (r) ( ir+ / J(t)N(x, t; s(r), T)g{r)dr = 0, 

J o 

where gr(i) = u(s(t),t), < t < T. Taking into account that iV^ = — G x and 

passing to a limit as x — > s(t) — 0, we obtain by Lemma 14.11 

(4.24) 

f G x (s(t),t;s(r),T)g(T)dT+ [ s'(r)N( S (t),t;s(r),T)g(T)dT = 0. 



We are going to explain that this integral equation for g has only the trivial 
solution g(t) = 0. One can easily see that for any T\ < T there is a constant 
C > such that for r G [0, T±] 

\G x (s(t),t-s(r),T)\ < -==, | S '(r)iV( S (t),t; S (r),r)| < ( 



l „ v .^ v „ VVJ „,„ v .„. /l _^__ 

Now, by (glgj) it follows that 

sup \g(t)\ < 8CVhsup \g(t)\, t x G (0,Ti). 
[0,ti] [0,ti] 
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Choose t\ so that %C^ft{ < 1; then we have g(t) = for < t < t\. 

The same argument shows that if g(t) = for t G [0,ioL then there is a 
5 > such that g(t) = for t E [0, t + 5]. Hence, g(t) = for t G [0, T), i.e., 
(|3.5p holds. This completes the proof of Lemma 14.21 □ 

5. Local existence-uniqueness 

We study the local existence-uniqueness properties of the system of non- 
linear integral equations (14. lip . (I4.12p by employing the Banach Contraction 
Fixed Point Theorem. 

Let s = const > 0, and let E be the space of all pairs of continuous 
functions (v(t),s'(t)), t G [0, e\. Equipped with the norm 

|| s 7 ) || £ = max{||w|| ej ||s'|| £ }, where ||v|| e = sup \v(i)\, ||s'|| £ = sup \s'(t)\, 

[0,e] [0,e] 

E is a Banach space. 

We fix a constant T > and introduce the norms 

\\f\\ T = sup \f(t)\, II/Ht = sup \f(t)l \\<p% = sup \$'(x)\. 

[0,T] [0,T] [0,6] 

In the following we may assume that e < T. 
Consider in E the operator 

(5.1) <S>(v, S ') = (A(v,s'),B(v,s')), 
where 

(5.2) A(v,s')(t) = 2A 1 (v,s')(t) + 2A 2 (v,s')(t) + 2A 3 (v,s')(t) for t > 
with 

A 1 (v,s')= [ N x {s{t),t;s(T),r)v(T)dT, A 2 (v,s') = -[ N x (s(t),t;0,r)f(r)dT, 



o 



Jo 



and 

r s(t) 

(5.3) B(v,s')(t)= W(v,s')(x,t)dx, 

Jo 

with 
(5.4) 



W(v, s')(x, t) = f(t)e~ xt - a + e~ xt V [* J v {v, t)d£, 

u=l Jo 



where J u are the integrals introduced in (14. 7p . 

In the above notations the system of integral equations (14. lip and (|4,12p 
could be written as 

(5.5) (v,s) = $(u,s'). 

Next we prove that locally (say, for < t < e) the equation (I5.5P has a unique 
solution by showing that for every large enough M > there is an e > such 
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that the operator maps the closed ball Bm = {(v,s') : \\(v,s')\\ £ < M} 
into itself, and its restriction on 6m is a contraction mapping. 
We impose the following a priori conditions on M and e : 

(5.6) M>1, e < min{l,6/(2M)}. 

Then ||s'|| £ < M implies \s(t) - b\ < Me < 6/2. Therefore, 



(5.7) 



\s'\\ £ <M =>• 6/2 < s(t) < 36/2 for < t < e. 



Straightforward computations show that the operator $ maps Bm into itself 
if 



(5.8) 
and 

(5.9) 



M > 1 + (8 + 56 2 )||^|| fe + 46 2 ||/|| T + 46 2 a 



£ < 



4(M + 4/6) 2 + (32/6)' 



T + 56 2 M + 2||/|| T 



In order to prove that the operator <3? : Bm — > Bm is a contraction map- 
ping we investigate the contraction properties of integral operators in (|5.2p - 
(|5.4p . Since A\(v,s') and A 2 (v,s') are Voltera type integral operators, one 
can easily prove that there exists e > such that for t G [0, e] 



1 



(5.10) lA^s'Jit) -A v (y 2 ,s' 2 )(t)\ < -11(^1-^2,4 



1,2, 



whenever (vi,Si), (v 2 ,s 2 ) G Bm- 

Next we consider A 3 (v, s'). Since G(x, t; £, 0) = K(x, t; £, 0)— K{— x, t; £, 0), 
we have A^(vi, Si)(t) — A$(v 2 , s' 2 )(t) = Aj(t) — A^(t), where 

rb 



1 



(n(*)±«r 



_(£2W±Ol 

e 4t 



£'(£K- 



2^ 7o 

By the elementary inequality | e~ xi —e~ X2 | < \x\—x 2 \e~ mm ^ Xl,X2 \ x±,x 2 > 
0, the expression in the square brackets in does not exceed by absolute 
value 

l*l(*)-«2(*)l (ai(t) + 8a(t) + 20 exp f _ J_ min[(si(t) + ^ (s2(t) + ^2- 



56 

T 



<-e- fe2/16i -|| S ' 1 -4l| £ 



because 6/2 < Si(t) < 36/2, i = 1,2 (see (|5.7jl ). Therefore, taking into 
account that e -62 / 16 * < 16i/6 2 , we obtain 

\A+(t)\ < lOvWllbK - s 2 || £ < - s 2 || £ if t < e < (80||^|| b )- 2 . 

In order to estimate A^ we write it in the form A^(t) = A^^{t) + A^ 2 (t), 
where 



^3,1 (t) 



1 



Jo 



b—8 , ,,2 
(fiW— gj 



1 — exp 



{si{t)-i? {s 2 {t)-if 



4t 



it 
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1 — exp 



At 



At 



2\pwt Jb-8 
and 5 = 2Mt 1 / 4 . 

From (|5.7p it follows that \si(t) + S2(t) — 2£| < 36 which implies 

-= ^ |ai(t)+a a (t)-2f| < — llsi-s 



At 



At 



1 °2\\e- 



Therefore, by the inequality \e x — 1| < l^le' 11 ', the expression in the square 
brackets in A 3 ~ 1 (t) does not exceed by absolute value (36/4) exp(36M/2)||s' 1 — 
S2 Thus, it follows 

r-6-5 



< ^ e 36M/2||~/|| 6 .|| s / 



s 2lle 



1 



(n(*)-fr 







2\/i 



Performing the change of variable z = in the latter integral, and 

estimating from above the resulting integral, we obtain 

h-(S-s 1 (i) 



/_ 



2V~t b—S {b-5- sl (t)) 

e dz < — e 

2Vt 



2 h M 2 

< — =e~^ < 8b -t 1/2 



2y/i 2^t 

z v c 

because 

5 + (ai(t) - 6) > 2Mt 1/4 - ||si|| £ • t > 2Mt 1/4 - Mt > Mt 1/4 , 
and (by flO]) and ([51)]) ) t" 1 exp(-M 2 /4 v / t) < 16. Therefore, 



A 31 {t) < 66 2 exp(36M/2)||^'|| b • ||si 



5 2 lie 



■Vt. 



Next we estimate A 3 2 {t). If £ G [6 — 5, 6], then (since si(0) = £2(0) = b) 

\sxit) + s 2 (t) - 2£| = |ai(t) - si(0) + s 2 (t) - s 2 (0) + 2(6 - £)| 

< (114 lie + II4IU) -t + 25< 2Mt + AMt l/i < 6Mi 1/4 , 
which implies 

(si(t)-e) 2 (S2(t)-t) 2 



At 



At 



\si(t) - s 2 (t)\ 
At 



ai (t)+ 82 (t)-2S\ < -Mt^Ws^-s'^. 



Estimating the expression in the square brackets in A 3 2 {t) as in the case of 
A 3l {t) and taking into account (|4,3p . we obtain 

A^ 2 {t)\ < ^Mt 1 / 4 exp(3M 2 )||^|| fc || S / 1 -a 2 || e . 

Thus, 

I A3 (t)\ < [66 2 exp(36M/2) + 3Mexp(3M 2 )] W\\ b \\s[ - s 2 \\ £ ■ t 1/4 , 
which implies that 

1^3 (*)l < l\Wl-S 2 \U 

if t < e < 8~ 4 [66 2 exp(36M/2) + 3Mexp(3M 2 )]~ 4 ||^'||^ 4 . 
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Now, the estimates for A 3 (t) and At (t) imply that if 
(5.11) 

e < min |80||^'|| 6 )~ 2 , 8~ 4 [6b 2 exp(36M/2) + 3M exp(3M 2 )] 4 ||^'||^ 4 } , 
then 

(5.12) \A 3 ( Vl ,s[)(t) - A 3 (v 2 ,s' 2 )(t)\ <^\\s[ - s' 2 \\ £ . 

The estimates (15.10P hold if e satisfies inequalities similar to (15. 9h and 
(|5.1ip . Moreover, one can prove that the operator B is a contraction in Bm 
if e satisfies similar restrictions. We omit the details, but it is important to 
note that the right-hand sides of (|5.9p , (|5.1ip and the analogous inequalities 
(which guarantee that the operator $ is a contraction on Bm with a contrac- 
tion coefficient < 1) are given by expressions that decrease if the parameters 
involved (such as b, M, \\<p'\\b, ||/||r, ll/llr) increase. 

Therefore, applying the Banach Contraction Fixed Point Theorem, we 
obtain the following statement. 

Lemma 5.1. (a) For each constant M > 1 which satisfies \5. 8\) there is a 
constant e > such that the system of integral equations j^. 1 i[ ) and $4-12 ) 
(with s(t) = b + §ls'(r)dT) has a unique solution (v(t), s'(t)), t 6 [0,e], 
such that \\v\\ e < M and \\s'\\ £ < M. 

(b) The constant e may be chosen so that 

(5.13) £ = ^,6,1/6,11/1^,11/11^11^116), 

where h(yi,y 2 ,y3,y4,y5,y6), y% > 0, is a monotone decreasing function with 
respect to each argument yi, i = 1, . . . , 6. 

Next we prove uniqueness of solutions of Problem P. 

Lemma 5.2. For each T < oo, Problem P has at most one solution for 
i€[0,T). 

Proof. Suppose that (ui(x, t), si(t)) and (u 2 (x,t), s 2 (t)) are two solutions of 
Problem P on the interval [0, T), T < oo. Then, in view of Lemma l4.2( the 
pairs of functions (v%(t), s'i(t)) and (v 2 (t), s' 2 (t)), where t>i(t) = ^-(si(t),t), 
v 2 (t) = ^jr(s 2 (t), t), are solutions of the system of integral equations (14. lip . 
(f37T2~|) . Fix e <T and choose M > 1 so that 

M > max{||ui(i)|| eo , K|| £o , ||u2(t)|| 6o , KlU}- 

By Lemma 15.11 there is a positive constant e < £q such that the pairs 
(vi(t), Si(t)) and (v 2 (t) , s 2 (t)) coincide on the interval [0,e]. Therefore, the 
integral representation (|4.10p implies 

ui(x, t) = u 2 (x, t), si(t) = s 2 (t) if 0<t<£, 0<x< s(t). 

Having proved uniqueness for a small time interval < t < e, we can 
proceed in a similar way, step by step, to get uniqueness for all t > 0. Let 
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to < T be a positive number such that 

(5.14) ui(x,t) = u 2 (x,t), si(t) = s 2 (t) if < t < t , < x < s(t). 

Then (e~ A *°tti(x, t + t ), si(t + to)) and (e~ xt °u 2 (x,t + t ), s 2 (t + t )) are 
two solutions of Problem P on the interval [0,T — to), if considered with 
ft(t) = f(t + t ), 

hit) = e- xto f{t+t ), <pi{x) = e- A *°ni(x,t ) = e- xto u 2 {x,to), h = Sl (t ) = s 2 

instead of f(t), f{t),(p(x) and b. By the above argument, there is a constant 
£\ > such that 

u\{x,t) = u 2 {x,t), s\{t) = s 2 (t) if < t < to + ei, < x < s(t). 

Therefore, (|5.14p holds for each to < T, i.e. the solutions {u\{x, t), s\(t)) 
and (u 2 (x, t), s 2 (t)) coincide on [0,T). □ 

6. Existence of global solution 

Lemma 14.21 and Lemma 15 . 1 1 guarantee that the Problem P has a solution 
for < t < e for sufficiently small e > 0. In order to prove the existence of 
a global solution we need a priori estimates for s(t) and u x (x,t). 

By Lemma [2.21 there are constants C\ = C\{T) and C 2 = C 2 (T) such 
that 

(6.1) \s'(t)\ < Ci, l/C 2 < \s(t)\ <C 2 , < t < T. 

Lemma 6.1. Suppose that the pair of functions (u(x,t), s(t)) is a solution 
of Problem P for t G [0,T), < T < oo. Then 

(6.2) := sup \u x (x, t)\ < oo, 

D T 

where D T = {(x,t) : < x < s{t), < t < T}. 
Proof. It is enough to prove that 

(6.3) m := sup \v(t)\ < oo, 

[0,T) 

where v(t) = u x (s(t),t). Indeed, by (|4.15|) . u x {x,t) = Ylu=i^( x ^)^ wnere 
the integrals I u (x,t) are given by (I4.16p . 

We have I\{x,t) = I\^{x,t) + I\ j2 {x, t), where Ii t i(x,t) and I\ j2 (x,t) are 
given in (I4.19j) and (|4.20p . First we estimate |Ji ; i(a;,t)| : 

■ , . | / 1 /"* \s(t) - s(t)\ - (-w-^) 2 , 

<™ - (— / o '4 (f _ T) 3/2 e 4(t " T) ^ 

air 



— f= / / NO n e 4 (*— ) exp 
^io 4(t-r)3/ 2 



■ ( S ( t )_ x )2_ (s(r) _ x) 2 n 
4(t-r) 

From (|6.ip it follows that |s(t) — s(r)|/(t — r) < C\, so the first integral 
in the brackets does not exceed J C\/{2sJt - r)dr < C\\fT . By (|6.ip . the 
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expression in the square brackets in the integrand of the second integral can 
be estimated from above by 

lS f { ~_ S ^ )l \s(t) + s(r)-2x\<C 1 C 2 . 

Therefore, in view of (|4.2p . the second integral does not exceed e ClC2 . Hence, 
\h,i(x,t)\<m- (dVT+e ^ 2 ) . 



By (|6.ip we have C 2 1 < x + s(t) < 2C 2 for x < s(t). From these inequal- 
ities and (|4.2p it follows 

\h,2(x,t)\ < m- 2C| / K x {C^ l ,t-{),T)dT <m-2C$. 
Jo 

On the other hand (|4.2p and (|4.3p imply 

(6.4) |/ 2 (x,t)|<||/||T, |/ 3 (x,t)|<2||^|| 6 . 

Hence, 

sup | ^.(x,*) | <m - (c 1 ^ + e ClC2 +2C 2 2 ) + ||/|| T + 2||^|| 6 , 

i.e., ([631) implies ffO]) . 

Next we prove ljO]l . By (|iH]> . = 2 Y%=\ L{s{t), t), where the inte- 
grals I v (x,t) are given by (|4.16p . 

First we consider Ii(s(t), t) = J ' N x (s(t),t; s(t),t)v(t)cLt. Since iV(x, t; £, r) 
if (x, t; £, r) + K(—x, t; £, r), we have 



|AT/ ,. , , ,. \s(t)-s(T)\ _ (»(0-»(r)) 2 |a(t)+s(r)| _ i£W+£W)g 
l^(Kt),t;Kr),r)|<LU__^ e ^ + ^-A^ e ^ . 

From (|6.ip it follows that |s(i) — s(r)|/(t — r) < Ci, so the first term on 
the right in the above inequality is less than Ci/(4\/i — r). On the other 
hand, ([S3]) implies 2/C 2 + s(r) < 2C 2 . Therefore, in view of (IOD . 



4(t - r) 3 / 2 ~ 4(t - r) 3 / 2 ~ 2 2^t^ ' 

Thus, we obtain 

(6.5) \N x ( S (t),t;s(r),r)\ < (d/2 + g|) i , < r < t < T. 

zyr — r 

Let 5 6 (0,T) (later we will choose 5 sufficiently small), and let 
fi(t) = sup{|u(r)| : < t < t}. 
For te(T-6,T), ([63]) implies 

|/i(*(t),t)| < Im(T-S) [ T & \N x (s(t),t;s{T),T)\dT+fi(t) f \N x (s{t),t;s(r), 
Jo Jt-s 

< fi(T - 5) ■ (d/2 + Cf ) VT + fi(t) ■ (d/2 + Cl)VS. 
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Therefore, in view of (|4.1ip and (16. 4p . we obtain 

(6.6) \v(t)\ < iu(T-5)iC 1 + 2Cl)Vf+^t)iC 1 +2Cl)V6+2\\f\\ T +4\\ l p'\\ b 

for t G (T-5,T). 
Choose 5 so that 

(Cx + 2Cl)VS < 1/2. 

Then (j6.6j) implies 

fi(t) <2fi(T -5)(C 1 + 2Cl)Vf + 4\\f\\ T + 8\\^\\ b for t£(T-5,T). 

Hence, m = sup{/i(i), i G [0, T)} < oo, i.e., (|6,3p holds, which completes 
the proof of Lemma 16.11 □ 

Proof of Theorem VS ."A By Lemma l5.2l Problem P has at most one global 
solution. Now we prove that Problem P has a global solution. Assume the 
contrary, and let T be the greatest positive number such that Problem P 
has a solution for t G [0, T). Let the pair of functions u(x,t) and s(t) be a 
solution of Problem P for t G [0, T). 

For each to < T we can consider a "modified" Problem P with data 
(6.7) 

fl(t) = f(t + t), fi(t) = e- xto f(t + t), ^) = e- Xto u(x,t ), h = s(t ) 

instead of f(t),f(i),<p(x) and 6. By the local existence-uniqueness result 
given in Lemma 15. 11 for each M\ > 1 which satisfies 

(6.8) Mi > l + (8 + 56?)||^i|| 6l +46?||/i|| T + 4&fa 
and each e > with 

(6.9) e < fc(Mi,6i, l/6i, H/iIIt, ||/i||t, ll^'illbj, 

the "modified" Problem P has a solution (ui(aj, t), si(t)) for < t < e. 
Then, the pair (U (x , t) , S (t)) with 

m lt)= fet) 0<t<t S(t) = l S{t) °<*<*0 

|e At °ui(x,t-t ) to<i<^o + e' \ai(t-*o) t <t<t + e 

is a solution of Problem P for < t < to + e. 

Moreover, in view of the a priori estimates given in Lemma 12.21 and 
Lemma 16. 1\ by Lemma 15.11 we can choose one and the same e for every 
to < T. Indeed, let us set 

M = l+(8+5C7 2 2 )^T+4C 2 2 ||/|| 2T +4C7 2 2 a, g = h(M, C 2 , C 2 , \\f\\ 2 T, Wfhr, *t). 

Therefore, choosing to > T — e, we get the existence of a solution of 
Problem P for t G [0, to + e) with to + i > T, which contradicts the choice of 
T. Hence Problem P has a global solution for t G [0, oo), i.e., Theorem 13.21 
holds. 

In view of Lemma 13. 1 1. this implies that Problem P has a unique global 
solution, i.e., Theorem 11.11 holds as well. 
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7. Concluding remarks 

1. During the last 40 years various mathematical models for evolution 
of tumors have been developed and analyzed - see the survey papers [U E] 
and the bibliography therein. Some of those models are in the form of free 
boundary problems for partial differential equations, whereby the tumor 
surface is a free boundary and the tumor growth is determined by the level 
of a diffusing nutrient concentration |1Q|. Illj (see also pQ— [S], [I])- The main 
physical and biological concepts underlying such type of models are the mass 
conservation law and reaction-diffusion processes within the tumor. Usually 
additional geometric assumptions on the shape of the tumor are imposed - 
see, for instance, [TOj, [8], where the tumor is supposed to be spherically 
symmetric. 

A slight modification of those models is considered in [T5]. It describes 
the growth of an avascular solid tumor which receives nutrient supply via 
a diffusion process only through some part of its boundary (called base of 
the tumor), and it is assumed that there is no nutrient flow through the 
remaining part of the boundary. Moreover, the tumor is supposed to be 
thin and approximately disc-shaped, so only one spatial dimension, say x, is 
considered. With tumor's base situated at x = the nutrient concentration 
a(x,t) satisfies the reaction-diffusion equation 

n da d 2 a , . . 

(7.1) c— = -^3- Xa, 0<x<s{t), t>0, 

where s(t) > is the tumor's thickness at time t, A = const > 0, Act is the 
nutrient consumption rate, and c > is a dimensionless constant coming as 
a ratio of the nutrient diffusion time scale to the tumor growth time scale. 

Following [10], it is assumed that all tumor cells are physically identical 
in volume and mass, and that the cell density is constant throughout the 
tumor. As in [8], the cell proliferation rate within the tumor is given by 
P(a) = fi(a — <?), where \x and a are positive constants. These assumptions 
lead to the equation 

(7.2) s '(t) = u {a(x,t) - a)dx, t > 0. 

J o 

In addition, the following mixed type boundary conditions hold: 

(7.3) a(Q,t) = f(t), t>0, 

(7.4) a(x, 0) = <p(x), x G [0, 6], s(0) = b > 0, y>(0) = /(0), 

(7.5) -£(s(t),t) = 0, t>0, 

where f{t) > is the external nutrient concentration at the base of the 
tumor at time t, <p(x) > is the initial nutrient concentration within the 
tumor, and the condition (|7.5h comes because it is assumed that there is no 
nutrient transfer through the free boundary x = s(t). 
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The parameters c, \i and A in (j7.1j) - (|7.5p depend on the choice of time 
and length units. One may scale out x and t in an appropriate way in order 
to get c = 1 and u = 1 (A may change as well), which shows that Problem 
(|7.ip ~ (j7.5p is equivalent to Problem P. 



2. Another interesting question in the study of mathematical models 
of tumor growth is under what conditions does the tumor grow, shrink or 
become dormant. In order to answer that question one needs to find the sta- 
tionary solution (which gives the dormant case) and analyze its asymptotic 
stability (see [H [6] and the bibliography there) . 

In the case of Problem P, if f(t) = a = const then it is easy to see that 
the stationary solution is given by the pair (u(x),b), where 



cosh ( v X(x — b) 

(7.6) u(x) = a 



cosh yy/Xb 
and b is determined by the equation 
(7.7) a tanh (bVX) =abV\. 



3. In (|1.6p of Theorem ll.lt the assumptions f(t) > for t G [0, oo) and 
f(x) > for x G [0, b] come from the corresponding mathematical model 
(Section 7.1). However, the result stated in Theorem 11.11 remains valid 
without those requirements. 

Indeed, let f(t) G C 1 ([0, oo)) and ip(x) G C 2 ([0, b]) be arbitrary functions. 
Then, under the assumptions of Lemma [2.21 the following a priori estimates 
hold: 

(7.8) \u(x,t)\<C T , 0<x<s{t), 0<t<T, 

(7.9) 

- (C T + a)s(t) < s'(t) < (C T - a)s(t), be~ {CT+li)t < s(t) < be {CT -~ a) \ 

where b = s(0), Ct = max ^sup[ 0T ) |/(t)|, sup[ 0j &] • The proof of 

Theorem II .11 is the same, but one needs to use the estimates (|7.8p and (|7.9p 
instead of (JOJ and (J23D in Lemma 

4. It is known that the free boundary in the one-dimensional Stefan 
problem (see Ch. 8], 02], [21 Ch. 17]) is a C°°-curve (see O H E3] and 
the bibliography therein). In the context of tumor models a similar result 
is proven in [6l Theorem 4.1]. 

In the case of our Problem P it is easy to see that s(t) G C 2 ([0,oo)). 
Indeed, since ut(x,t) is defined and continuous for < x < s(t), t > 0, from 
(Ojl it follows 

s"(t) = [u(s{t),t)-a]s'{t)+ / u t (x,t)dx. 

Jo 
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Therefore, taking into account that 




s"(t) = [u(s(t),t) - a- \]s'(t) - u x (0,t) - Xas(t), 



where the expression on the right is a continuous function for t > 0, i.e., 



However, higher derivatives of s(t) may not exist if we assume / G 
C 1 ([0, oo)) only, since the condition fjl .21) in Problem P is nonlocal (compare 
with the case of one-dimensional Stefan problem, where the infinite differ- 
entiability of the free boundary does not require infinite differentiability of 
the boundary data at x = - see [E]). In our case, one can prove the fol- 
lowing: In Problem P, the free boundary x = s(t), t £ (0, oo) is an infinitely 
differentiable curve if and only if f(t) G C°°((0,oo)). We will present the 
details somewhere else. 
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